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Abstract. Traditional plasma physics has mainly focused on regimes 
characterized by high temperatures and low densities, for which quantum- 
mechanical effects have virtually no impact. However, recent technolog- 
ical advances (particularly on miniaturized semiconductor devices and 
nanoscale objects) have made it possible to envisage practical applica- 
tions of plasma physics where the quantum nature of the particles plays 
a crucial role. Here, I shall review different approaches to the modeling 
of quantum effects in electrostatic collisionless plasmas. The full kinetic 
model is provided by the Wigner equation, which is the quantum analog 
of the Vlasov equation. The Wigner formalism is particularly attractive, 
as it recasts quantum mechanics in the familiar classical phase space, 
although this comes at the cost of dealing with negative distribution 
functions. Equivalently, the Wigner model can be expressed in terms 
of N one-particle Schrodinger equations, coupled by Poisson's equation: 
this is the Hartree formalism, which is related to the 'multi-stream' ap- 
proach of classical plasma physics. In order to reduce the complexity of 
the above approaches, it is possible to develop a quantum fluid model by 
taking velocity-space moments of the Wigner equation. Finally, certain 
regimes at large excitation energies can be described by semiclassical 
kinetic models ( Vlasov- Poisson), provided that the initial ground-state 
equilibrium is treated quantum-mechanically. The above models are 
validated and compared both in the linear and nonlinear regimes. 



1 Introduction 

Plasma physics deals with the TV-body dynamics of a system of charged par- 
ticles interacting via electromagnetic forces. The study of plasmas arose in the 
early twentieth century when scientists got interested in the physics of gas dis- 
charges. After World War II, plasmas became the object of intensive experimental 
and theoretical research, mainly because of the potential applications of nuclear 
fusion, both military (hydrogen bomb) and peaceful (energy production through 
controlled thermonuclear fusion). In parallel, plasma physics was also developed 
by astrophysicists and geophysicists, which is not surprising, as it is thought that 
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about 90% of all matter in the visible universe exists in the form of a plasma. More 
precisely, plasmas are observed in the Sun surface, the earth's magnetosphere, and 
the interplanetary and interstellar media. 

Both fusion and space plasmas are characterized by regimes of high tempera- 
ture and low density, for which quantum effects are totally negligible . However, 
physical systems where both plasma and quantum effects coexist do occur in nature, 
the most obvious example being the electron gas in an ordinary metal. In metals, 
valence electrons are not attached to a particular nucleus, but rather behave as free 
particles, which is why metals make good electric conductors. Although some level 
of understanding of metallic properties is achieved by considering noninteracting 
electrons, a more accurate description can be obtained by treating the electron pop- 
ulation as a plasma, globally neutralized by the lattice ions. At room temperature 
and standard metallic densities, quantum effects can no longer be ignored, so that 
the electron gas constitutes a true quantum plasma. 

For ordinary metals, however, the properties of the electron population (band 
structure, thermodynamic properties) are mainly determined by the presence of a 
regular ion lattice, typical plasma effects being only a higher order correction. In 
recent years, though, there has been tremendous progress in the manipulation of 
metallic nanostructures (metal clusters, nanoparticles, thin metal films) constituted 
of a small number of atoms (typically 10— 10 5 ) [T ], l2l l3l 14), 13] . For such objects, no 
underlying ionic lattice exists, so that the dynamics of the electron population is 
principally governed by plasma effects, at least for large enough systems. Further, 
the development of ultrafast (femtosecond and, more recently, attosecond) laser 
sources makes it possible to probe the electron dynamics in metallic nanostructures 
on the typical time scale of plasma phenomena, which is indeed of the order of the 
femtosecond. Metallic nanostructures thus constitute an ideal arena to study the 
dynamical properties of quantum plasmas. 

Another possible application of quantum plasmas arises from semiconductor 
physics [fU HI 13 llOj - Even though the electron density in semiconductors is 
much lower than in metals, the great degree of miniaturization of today's electronic 
components is such that the de Broglie wavelength of the charge carriers can be 
comparable to the spatial variation of the doping profiles. Hence, typical quantum 
mechanical effects, such as tunneling, are expected to play a central role in the 
behavior of electronic components to be constructed in the next years. 

Finally, quantum plasmas also occur in some astrophysical objects under ex- 
treme conditions of temperature and density, such as white dwarf stars, where the 
density is some ten order of magnitudes larger than that of ordinary solids 
Because of such large densities, a white dwarf can be as hot as a fusion plasma (10 8 
K), but still behave quantum- mechanically. 

A classical system of charged particles qualifies as a plasma if it is quasineutral 
and if collective effects play a significant role in the dynamics |12j . 'Quasineutral- 
ity' means that charge separation can only exist on a short distance, which, for 
classical plasmas, is given by the Debye length. On distances larger than the Debye 
length, the plasma is basically neutral, except for small fluctuations. By 'collective 
effects' we mean particle motions that depend not only on local conditions, but 
also - indeed principally - on the positions and velocities of all other particles in 

1 Quantum effects do play an important role, for example, in determining the fusion cross 
sections. But this is a nuclear physics rather than a plasma physics issue. What we mean is 
that the dynamics and thermodynamics of fusion and space plasmas are completely unaffected by 
quantum corrections. 
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the plasma. (In solid state and nuclear physics, collective effects are usually called 
'mean-field' effects, because they arise from the average field created by all the 
particles). Such collective behavior is possible because of the long-range nature 
of electromagnetic forces. In contrast, for neutral gases, the dominant interaction 
mechanism is provided by short-range molecular forces of the Lennard- Jones type: 
individual molecules thus move undisturbed in the gas until they make a collision 
with another molecule (which occurs when the interaction potentials overlap). It 
should be added that real plasmas are often only partially ionized, so that a frac- 
tion of neutral molecules is also present. In order to have true plasma behavior, 
we must therefore require that the collision rates of electrons and ions with the 
neutral molecules be relatively low compared to typical collective phenomena. In 
the present work, we shall avoid this issue altogether by restricting our analysis to 
the simpler case of fully ionized gases. 

When quantum effects start playing a role, the above picture gets more compli- 
cated, as an additional length scale is introduced, namely the de Broglic wavelength 
of the charged particles, Xb = h/mv. The dc Broglie wavelength roughly repre- 
sents the spatial extension of the particle wave function - the larger it is, the more 
important quantum effects are. From the definition of Xb, it is clear that quan- 
tum behavior will be reached much more easily for the electrons than for the ions, 
due to the large mass difference. Indeed, in all practical situations, even the most 
extreme, the ion dynamics is always classical, and only the electrons need to be 
treated quantum-mechanically. In the present paper, we shall always refer implic- 
itly to electrons when discussing quantum effects. In addition, only electrostatic 
(Coulomb) interactions will be considered. Magnetic fields do introduce novel and 
interesting effects, but the fundamental properties of quantum plasmas are already 
present in the purely electrostatic scenario. 

In the rest of this paper, we shall first obtain a number of qualitative results by 
using simple arguments from dimensional analysis. This will be useful to extract 
the relevant dimensionless parameters that determine the various physical regimes 
(classical/quantum, collisionless/collisional). Subsequently, we shall derive and il- 
lustrate several mathematical models that are appropriate to describe the dynamics 
of a quantum plasma in the collisionless regime. 



2 Physical regimes for classical and quantum plasmas 

In this section, we shall derive a number of parameters that represent the typi- 
cal length, time, and velocity scales in a classical or quantum plasma. These can be 
obtained using elementary considerations based on dimensional analysis. Of course, 
more detailed studies would be necessary to understand how such parameters ac- 
tually intervene in real physical phenomena (for instance, whether a certain time 
scale represents a typical oscillation frequency, or rather a damping rate). Here, we 
shall derive the algebraic expression for these quantities and simply state, without 
proof, what they represent physically. 

In addition, it will also be important to establish certain dimensionless param- 
eters. Dimensionless parameters allow us to discriminate between different physical 
regimes, characterized by situations where one effect dominates over another. In 
particular, we shall look for parameters that define whether a plasma is classical or 
quantum, and whether it is dominated by individual effects (collisionat) or collective 
effects (collisionless). 
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2.1 Classical plasmas. We consider a plasma of number density n, composed 
of particles (typically, electrons) with mass m and electric charge e, interacting via 
Coulomb forces (hence the electric permittivity eq). With these four parameters, 
we are able to construct a quantity that has the dimensions of an inverse time, i.e. 
a frequency: 

1/2 

) ■ (2-1) 

me J 

The latter quantity is known as the plasma frequency and it represents the 
typical oscillation frequency for electrons immersed in a neutralizing background 
of positive ions, which is supposed to be motionless because of the large ion mass. 
The oscillations arise from the fact that, when a portion of the plasma is depleted 
of some electrons (thus creating a net positive charge) , the resulting Coulomb force 
tends to pull back the electrons towards the excess positive charge. Due to their 
inertia, the electrons will not simply replenish the positive region, but travel further 
away thus re-creating an excess positive charge. In the absence of collisions, this 
effect gives rise to undamped electron oscillations at the plasma frequency. 

Note that the plasma frequency is independent on the temperature. If we do 
introduce a finite temperature T, then we can construct a typical velocity: 

where ks is Boltzmann's constant. This is the thermal velocity, which represents, 
just like in ordinary gases, the typical speed due to random thermal motion. 

By combining the above two quantities, one can define a typical length scale, 
the Debye length: 

v T ( SoksT^ 



X D = -=r-^ ■ (2.3) 



The Debye length describes the important phenomenon of electrostatic screening: if 
an excess positive charge is introduced in the plasma, it will be rapidly surrounded 
by a cloud of electrons (which are more mobile and thus react quickly). As a result, 
the positive charge will be partially screened and will be virtually 'invisible' to 
other particles situated at a large enough distance. Quantitatively, this amounts to 
saying that the electrostatic potential generated by an excess charge does not fall, 
like in vacuum, as 1/r, but rather obeys a Yukawa-like potential exp(—r/\£>)/r, 
which of course decays much more quickly and on a distance of the order of the 
Debye length. The Debye screening is at the origin of one of the most crucial of all 
plasma properties, namely quasineutrality: charge separation in a plasma can exist 
only on scales smaller than Ad, but it is screened out at larger sales. 

Let us now try to construct a dimcnsionless parameter using the above quanti- 
ties: m, e, £oi n, and T. It is easily seen that only one such parameter exists, and 
it reads as 

e 2 1/3 

9C = ■ (2.4) 

This is known as the (classical) graininess parameter or coupling parameter. 
It is illuminating to show that gc can be written as the ratio of the interaction 
(electric) energy E lnt to the average kinetic energy E^ n . Indeed, for particles 
situated at typical interparticle distance d = n^ 1 ^ 3 , one has E int — e 2 /(eod) and 
Skin = k B T, which immediately yields the expression (|2.4ll . 
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The expression gc = E- lnt /-Ekm allows us to guess the physical relevance of the 
coupling parameter. When gc is small, the plasma is dominated by thermal ef- 
fects, whereas two-body Coulomb interactions (i.e. binary collisions) remain weak. 
In this regime, the main held acting on the charged particles is the nonlocal mean 
field, which is responsible for typical collective effects. This is known as the colli- 
sionless regime. On the contrary, when gc — 1 or larger, binary collisions cannot 
be neglected and the plasma is said to be collisional or strongly coupled. We also 
note that, following lj2.4|l . classical plasmas are collisionless at high temperatures 
and low densities. 

Alternatively, gc can be written as the inverse of the number of particles con- 
tained in a volume of linear dimension \p>, raised to a certain power: 

/ 1 X2/3 

This shows that a plasma is collisionless when the Debye screening is effective, i.e. 
when a large number of electrons are available in a Debye volume. 

2.2 Quantum plasmas. Quantum effects can be measured by the thermal 
de Broglie wavelength of the particles composing the plasma 

A B = — , (2.6) 
mvT 

which roughly represents the spatial extension of a particle's wave function due to 
quantum uncertainty. For classical regimes, the de Broglie wavelength is so small 
that particles can be considered as pointlike (except, as mentioned in the Introduc- 
tion, when computing collision cross-sections), therefore there is no overlapping of 
the wave functions and no quantum interference. On this basis, it is reasonable to 
postulate that quantum effects start playing a significant role when the de Broglie 
wavelength is similar to or larger than the average interparticle distance n -1 / 3 , i.e. 
when 

n\% > 1 . (2.7) 

On the other hand, it is well known from the statistical mechanics of ordinary 
gases |13| that quantum effects become important when the temperature is lower 
than the so-called Fermi temperature Tp, defined as 

k B T F = E F = ^- (3vr 2 ) 2 / 3 n 2 / 3 , (2.8) 

where we have also defined the Fermi energy Ep. When T approaches Tp, the 
relevant statistical distribution changes from Maxwell-Boltzmann to Fermi-Dirac. 
Now, it is easy to see that the ratio \ = T /Tp is simply related to the dimensionlcss 
parameter n\^ B discussed above: 

X^^i^ 2 ) 2 / 3 ^) 2 / 3 . (2.9) 

Thus, quantum effects become important when x > 1- 

We now want to establish the typical space, time, and velocity scales for a 
quantum plasma, as well as the relevant dimensionless parameters. First of all, 
we stress that simple expressions can be found only in the limiting cases T ^> Tp 
(corresponding to the classical case treated previously) and T -C Tp, which is 
the 'deeply quantum' (fully degenerate) regime that we are going to analyze. Of 
course, there will be a smooth transition between the two regimes, but this cannot 
be treated using straightforward dimensional arguments. 
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Concerning the typical time scale for collective phenomena, this is still given by 
the inverse of the plasma frequency (|2.1(l . even in the quantum regime. However, the 
thermal speed becomes meaningless in the very low temperature limit, and should 
be replaced by the typical velocity characterizing a Fermi-Dirac distribution. This 
is the Fermi velocity: 

vp = (^M 1/2 = 1 (3 ^2 n) l/3 (210) 



With the plasma frequency and the Fermi velocity, we can define a typical 
length scale 

Af = — i (2.11) 

which is the quantum analog of the Debye length. Just like the Debye length, Af 
describes the scale length of electrostatic screening in a quantum plasma. 

The quantum coupling parameter can be defined as the ratio of the interaction 
energy E- mt to the average kinetic energy Eun- The interaction energy is the same 
as in the classical case, whereas the kinetic energy is now given by the Fermi 
energy Ey_- m = Ep. With these assumptions, one can write the quantum coupling 
parameter as 

E- mt 2 e 2 m / 1 \ 2 ^ 3 ( ftw p ^ ' 



9Q ~ E F (3tt 2 ) 2 /3 tPe nV3 ~ [ n \ F ) ~\E F ) ' (2 ' 12) 

where we have left out proportionality constants for sake of clarity. The third 
expression of gg in (|2.12|) is completely analogous to the classical one when one 
substitutes Af — > A rj . The last expression is more interesting, as it has no classical 
counterpart: it describes the coupling parameter as the ratio of the 'plasmon energy' 
hu) p (energy of an elementary excitation associated to an electron plasma wave) to 
the Fermi energy 2 . 

The quantum collisionless regime (where collective, mean-field effects dominate) 
is again defined as the regime where the quantum coupling parameter is small. From 
(|2.12() . it appears that a quantum plasma is 'more collective' at larger densities, in 
contrast to a classical plasma [see (|2.4|) ]. This may seem surprising, but can be 
easily understood by invoking Pauli's exclusion principle, according to which two 
fcrmions cannot occupy the same quantum state. In a fully degenerate fermion gas, 
all low-energy states are occupied: if we add one more particle to the gas, it will 
necessarily be in a high-energy state. Therefore, by increasing the gas density, we 
automatically increase its average kinetic energy, which, in virtue of 12. 1211 , reduces 
the value of gg.. 

2.3 Plasma regimes. We have so far defined three dimensionless parameters 
that determine whether the plasma is classical or quantum, and, in either case, 
whether it is collisional or collisionless: 

1. x = Tf/T : classical/quantum 

2. gc ■ collisional/collisionless (classical regime) 

3. <7q : collisional/collisionless (quantum regime) 

These parameters are functions of the temperature and density. In Fig. ^ wc 
plot on a log T — log n diagram the straight lines corresponding to \ — gc = gQ = 1, 
which delimitate the various plasma regimes |15j . 



2 gQ is directly proportional to the parameter r s /ao (where ao is Bohr's radius), commonly 
used in solid state physics 1141 . 
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Figure 1 Plasma diagram in the logT — logn plane. IONO: ionospheric 
plasma; SPACE: interstellar space; CORONA: solar corona; DISCHA: typi- 
cal electric discharge; TOK: tokamak experiment (magnetic confinement fu- 
sion); ICF: inertial confinement fusion; MET: metals and metal clusters; JUP: 
Jupiter's core; DWARF: white dwarf star. 



The logT — logn plane is divided into four regions, two of which are classical 
(above the T = Tp line) and two quantum. Each quantum/classical region is 
then divided in two collisional/collisionlcss subrcgions, identified by the kinetic 
equations that are relevant to each regime. As we shall see in the forthcoming 
sections, the Vlasov and Wigner equations are the appropriate models respectively 
for collisionless classical and collisionless quantum plasmas. 'Boltzmann' is used as 
a generic term for collisional kinetic equations in the classical regime. Collisional 
effects in the quantum regime are much harder to deal with, and no uncontroversial 
kinetic model exists in that regime, which is identified by 'Wigner (+ coll.)' on the 
figure. 

We point out that all previous considerations have implicitly assumed thermal 
equilibrium. Out-of-cquilibrium regimes should be treated much more carefully and 
the above results may not be entirely correct. For example, if an electron beam 
is injected into a plasma, the beam velocity will have to be taken into account 
when determining the de Broglie wavelength, which will therefore be smaller. For 
this reason, systems that are far from equilibrium can sometimes be treated with a 
semiclassical model, even though the corresponding equilibrium may still be fully 
quantum (see Sec. 14.411 . 

Several points corresponding to natural and laboratory plasmas have also been 
plotted in Fig. 1. We note that space and magnetic fusion plasmas fall in the clas- 
sical collisionless region, whereas inertial confinement fusion plasmas may display 
quantum and/or strong coupling effects. Extremely dense astrophysical objects 
such as white dwarf stars are definitely quantum and collisionless, even though 
they are as hot as fusion or solar plasmas. 
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Figure 2 Fcrmi-Dirac energy distribution for a case with T/Tp = 0.1. The 
vertical lines define a shell of thickness fcgT around the Fermi surface. 



3 Electrons in metals and metallic nanostructures: Pauli blocking 

The typical quantum coupling parameter for ordinary metals is larger than 
unity, so that, in principle, electron-electron collisions are as important as collec- 
tive effects. If that were really the case, one should abandon collisionless models 
altogether and resort to the full iV-body problem. This is hardly a feasible task. 
Fortunately, however, the effect known as Pauli blocking reduces the collision rate 
quite dramatically in most cases of interest. This occurs when the electron distri- 
bution is close to the Fermi-Dirac equilibrium at relatively low temperatures. The 
fundamental point is that, when all lower levels are occupied, the exclusion principle 
forbids a vast number of transitions that would otherwise be possible |14) . In par- 
ticular, at strictly zero temperature, all electrons have energies below Ep, and no 
transition is possible, simply because there are no available states for the electrons 
to occupy. At moderate temperatures, only electrons within a shell of thickness 
ksT about the Fermi surface (i.e. the region where E = Ep) can undergo collisions 
(this shell is delimited by the two vertical lines in Fig. For such electrons, the 
e-e collision rate (inverse of the lifetime r ee ) is proportional to fegT/fi (this is a 
form of the uncertainty principle, energy x time = const.). The average collision 
rate is obtained by multiplying the previous expression by the fraction of electrons 
present in the shell of thickness T about the Fermi surface, which is ~ T/Tp. One 
obtains 

k R T 2 , . 

(3.1) 



In normalized units, this expression reads as 



E F / T y _ 1 / T 



(3.2) 



Thus, u ee < lu p in the region where T < Tp and <?q > 1, which is the relevant 
one for metallic electrons (Fig. QJ. Restoring dimensional units, we find that, at 
room temperature, r ee ~ 10~ 10 s, which is much larger than the typical collisionless 
time scale t p = 27toj~ 1 ~ 10~ 15 s = 1 fs. In addition, the typical time scale for 
electron-lattice collisions T e i ~ 10~ 14 s = 10 fs is also larger than r p . Therefore, it 
appears that a collisionless regime is indeed relevant on a time scale of the order of 
the femtosecond. 
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A word of caution is in order, however, not to overestimate the effect of 
Pauli blocking. The above considerations are valid at thermodynamic equilibrium, 
whereas many more transitions are allowed for out-of-equilibrium electrons. There- 
fore, for strongly excited systems where many nonequilibrium electrons are present, 
the collision frequency may be larger than the simple estimate given in (|3.2I) . 

Some typical parameters for metallic electrons (gold) at room temperature are 
summarized in Table 1. We note that t p , the typical collisionless time scale, is of the 
order of the femtosecond. With the recent development of ultrafast laser sources 
with femtosecond period, it is therefore possible to probe the mean-field properties 
of metallic nanostructures. For example, the electron dynamics in thin gold films 
excited with femtosecond lasers was studied experimentally in several works [2115], 
We also point out that Xp, the typical collisionless length scale, is of the order of 
the Angstrom, which is comparable to the typical atomic size. 

Table 1 Typical parameters for electrons in gold at room temperature 



n 


5.9 x 10 28 m- 3 


T 


300 K 


LU p 


1.37 x 10 16 s- 1 


Tp 


0.46 fs 


Tie 


30 fs 


T F 


6.4 x 10 4 K 


E F 


5.53 eV 


V F 


1.4 x 10 b ms" 1 


Af 


1 x 10- 10 m 


9Q 


12.7 


r s /a 


3 



4 Models 

The most fundamental model for the quantum ./V-body problem is the Schrodin- 
ger equation for the ./V-particle wave function ip(xi,X2, ■ ■ ■ ,Xn,€). Obviously, this 
is an unrealistic task, both for analytical calculations and numerical simulations. 
A drastic, but useful and to some extent plausible, simplification can be achieved 
by neglecting two-body (and higher order) correlations. This amounts to assume 
that the ./V-body wave function can be factored into the product of N one-body 
functions: 

ip(x 1 ,x 2 , ...,x N ,t) = ^i(xx,t) ip2(x2,t) ■ ■ ■ ipir(x N ,t). (4.1) 

For fermions, a weak form of the exclusion principle is satisfied if none of the wave 
functions on the right-hand side of (|4.1() are identical 3 . 

The assumption of weak correlation between the particles is satisfied, as we 
have seen, when the quantum coupling parameter gQ is small. The set of N one- 
body wave functions is known as a quantum mixture (or quantum mixed state) and 

3 A stronger version of the exclusion principle requires that ip(xi, X2, . . . , £jv, t) is antisym- 
metric, i.e. that it changes sign when two of its arguments are interchanged. This can be achieved 
by taking, instead of the single product of N wave functions as in 14.11 , a linear combinations of 
all products obtained by permutations of the arguments, with weights ±1 (Slater determinant) 
1141 . This is at the basis of Fock's generalization of the Hartree model, as pointed out in Sec. 14.21 
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is usually represented by a density matrix 

N 

p(x,y,t) = ^2p a ip a {x,t) ip a {y,t), (4.2) 

where, for clarity, we have assumed the same normalization f\ip a fdx — 1 for all 
wave functions and then introduced the occupation probabilities p a . 

Both the Wigner and Hartree models described below are completely equivalent 
to models based on the density matrix formalism (Von Neumann equation). 

4.1 Wigner-Poisson. The Wigner representation |16| is a useful tool to ex- 
press quantum mechanics in a phase space formalism (for reviews see |17U18lll9| V 
As detailed above, a generic quantum mixed state can be described by N single- 
particle wave functions tjj a (x,t) each characterized by a probability p a satisfying 
^2 a= iPa — 1- The Wigner function is a function of the phase space variables (x, v) 
and time, which, in terms of the single-particle wave functions, reads as 



f(x, v, f ) = E ^ P° / + °° ^* (* + ^ ( x ~ e ™ UVri dX ( 4 - 3 ) 

(we restrict our discussion to one-dimensional cases, but all results can easily be 
generalized to three dimensions). It must be stressed that the Wigner function, 
although it possesses many useful properties, is not a true probability density, as 
it can take negative values. However, it can be used to compute averages just like 
in classical statistical mechanics. For example, the expectation value of a generic 
quantity A(x, v) is defined as: 

(A\ = I L f(x,v)A(x,v)dxdv 

[ ' JJf{x,v)dxdv ' { ' 

and yields the correct quantum-mechanical value . In addition, the Wigner func- 
tion reproduces the correct quantum-mechanical marginal distributions, such as the 
spatial density: 

n(x,t)= / f(x, v, t) dv = ^2 p a \ il> a | 2 . (4.5) 

a=l 

We also point out that, of course, not all functions of the phase space variables 
are genuine Wigner functions, as they cannot necessarily be written in the form 
(14.3(1 . In general, although it is trivial to find the Wigner function given the N wave 
functions that define the quantum mixture, the inverse operation is not generally 
feasible. Indeed, there are no simple rules to establish whether a given function of 
x and v is a genuine Wigner function. For a more detailed discussion on this issue, 
and some practical recipes to construct genuine Wigner functions, see |2U| . 

The Wigner function obeys the following evolution equation: 

9£ + v dl + 
dt dx 



2iirh 2 



dX dv'e im(v - v ' )x/h 



X ~2 



f(x,v',t) = 0, (4.6) 



4 For variables whose corresponding quantum operators do not commute (such as xv). 14.41 
must be supplemented by an ordering rule, known as Weyl's rule 1191 . 
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where 4>(x,t) is the self-consistent electrostatic potential. Developing the integral 
term in 1)4.6(1 up to order 0(h 2 ) we obtain 

2l +v 2l-±^^l = J^^^l + (h i ) (4 7) 

dt dx m dx dv 24m 3 9a; 3 dv 3 

The Vlasov equation is thus recovered in the formal semiclassical limit H — ► 0. We 

stress, however, that rigorous asymptotic results are much harder to obtain and 

generally involve weak convergence. 

The Wigner equation must be coupled to the Poisson's equation for the electric 

potential 

d 2 S e 



dx2 - eoKJ fdv-n ) , (4.8) 

where we have assumed that the ions form a motionless neutralizing background 
with uniform density rig (this is known as the 'jellium' model in solid state physics). 

The resulting Wigner-Poisson (WP) system has been extensively used in the 
study of quantum transport O I21| . Exact analytical results can be obtained 
by linearizing ((4.6(1 and ((4.8(1 around a spatially homogeneous equilibrium given by 
fo{v). By expressing the fluctuating quantities as a sum of plane waves exp(ikx — 
iu)t) with frequency u> and wave number ft, the dispersion relation can be written 
in the form s(k,ui) = 0, where the dielectric constant e reads, for the WP system, 

e WP (w, ft) = H f / — — dv . (4.9) 



risk J hk(u) — kv) 

With an appropriate change of integration variable, l(4.9l) can be written as 

bwp(w 1 *) = 1-^/ ? , I° {V l, 4/A 2 dv . (4.10) 

no J [u> — kvy — a k I'kwr 

From ((4.9|) , one can recover the Vlasov-Poisson (VP) dispersion relation by taking 
the semiclassical limit h — » 

e V p(-,fc) = l + 4/^ d - 
noft J u> — kv 

The integration in 1(4.9(1 and ((4.11|) should be performed along the Landau 
contour in the complex (Re v, Im v) plane, so that the singularity at v = u>/k 
is always left above the contour 12 . This prescription allows one to obtain the 
correct imaginary part of e(k, lu), which is at the basis of the phenomenon of Landau 
damping. 

Just like in the classical case (Vlasov-Poisson) , the dispersion relation ((4.9(1 can 
also support unstable solutions, i.e. solutions with a positive imaginary part of the 
frequency. These solutions grow exponentially until some nonlinear effect kicks in 
and leads to saturation of the instability. The stability property of the dispersion 
relation of the Wigner-Poisson system have been extensively studied in (221 1231 . 



4.2 Hartree. A completely equivalent approach to the WP system is obtained 
by making direct use of the N wave functions if) a (x, t). These obey N independent 
Schrodinger equations, coupled through Poisson's equation 

dij} a h 2 d 2 tp a 

lH ^t = -2^^- e ^' a=1 -- iV (412) 

S - ^(E^l^l 2 --) • (4.13) 
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This type of model was originally derived by Hartree in the context of atomic 
physics, with the aim of studying the self-consistent effect of atomic electrons on 
the Coulomb potential of the nucleus. Subsequently, Fock introduced a correction 
that accounts for the parity of the ./V-particle wave function for an ensemble of 
fermions (Hartree- Fock model), but this development will not be considered in this 
paper 5 . 

Instead, it is useful to think of the above Schrodinger-Poisson equations l|4.12[l - 
(|4.13|) as the quantum- mechanical analog of Dawson's multistream model |26| . 
Dawson supposed that the classical distribution function can be represented as 
a sum of N 'streams', each characterized by a probability p a , a density n a , and a 
velocity u a : 

N 

f(x,v,t) = ^2 

Pa Tla. 

(x,t) 5(v — u a (x,t)) , (4-14) 

a=l 

where S stands for the Dirac delta. The streams represent infinitely thin filaments 
in phase space. If / obeys the Vlasov equation, then the functions n a and u a each 
satisfy the following continuity and momentum conservation equations: 

^ + = 0, (4.15) 

Ot ox 

du a du a e d(f> 

-^rr + u a^r~ = — jr- , (4.16) 

Ot ox max 

coupled of course to Poisson's equation with the electron density given by n(x, t) — 
"Yli a p a n a (x 1 t). Note that the representation l|4.14[) presents some drawbacks, as the 
functions u a (x,t) can become multivalued during the time evolution. This means 
that the system (|4. 151) — (|4. lfc>|l will develop singularities, such as an infinite density 
at certain positions. When this happens, the fluid description H4.15fl - I|4.16|) ceases 
to be valid, although the phase space picture of the streams is still correct. 

This line of reasoning can be extended to the quantum case |27U28) by applying 
the Madelung representation of the wave function to the system (I4.12|) - (I4.13|) . Let 
us introduce the real amplitude A a (x,t) and the real phase S a (x,t) associated to 
the pure state ip a according to 

%jj a = A a exp(iS a /H) . (4.17) 

The density n a and the velocity u a of each stream are given by 

n a = \i> a \ 2 =A 2 a , Ua = -^. (4.18) 

m ox 

Introducing Eqs. (|4.17l) - (|4.18|) into Eqs. (|4.12|) - (|4.13l) and separating the real and 
imaginary parts of the equations, we find 

¥ + = °' (419) 

^a +u ^a _ e d<j> | h 2 d f d' 2 {^JWZ)/dx 2 \ 
dt a dx mdx 2m? dx \ ^/n^ J 

Quantum effects are contained in the ^-dependent term in (|4.20|) (sometimes called 
the Bohm potential). If we set h = 0, we obviously obtain the classical multistream 
model (I4.15|) (|4.1(j|) . An attractive feature of the quantum multistream model is 



5 A very similar model to the Hartree equations 14.121 — lA. 131 is known as TDDFT (time- 
dependent density functional theory) 1241 . Its linearized version goes under the name of RPA 
(random-phase approximation) 1251 . 
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that, contrarily to its classical counterpart, it does not generally develop singulari- 
ties. This is thanks to the Bohm potential, which, by introducing a certain amount 
of wave dispersion, prevents the density to build up indefinitely. 

Linearizing equations 14. f 9JI - 14.20J1 (supplemented by Poisson's equation) around 
the spatially homogeneous equilibrium: n a = no, u a = UQ a and 4> — 0, one obtains 
the following dielectric constant 

N LU 2 



a=l 



The classical multistream relation is obtained simply by setting h — in l|4.2I[l . 

The equivalence of the Wigner-Poisson and Hartree models can be readily 
proven on the linear dispersion relation. The homogeneous equilibrium described 
above corresponds to wave functions 

ipa = y/m exp fi mUaa x) . (4.22) 



h 

The Wigner transform (|4.3|1 of the above wave functions 1(4.221) is given by the 
following expression 

JV 

fo{v) = y^/p a n S(v - u 0a ), (4.23) 

a=l 

where S stands for the Dirac delta. By inserting ((4.23() into the WP dispersion 
relation l|4.9|l . one recovers precisely the multistream dispersion relation (|4.21|l . 

4.3 Fluid model. In classical plasma physics, fluid (or hydrodynamic) mod- 
els are often derived by taking moments of the appropriate kinetic equation (e.g. 
Vlasov's equation) in velocity space. The moment of order s is defined as: 

/ + OC 
f(x,v,t)v s dv. (4.24) 
-oo 

Then, the zeroth order moment is the spatial density and obeys the continuity 
equation; the first order moment is the average velocity and obeys a momentum 
conservation equation; the second order moment is related to the pressure, and so 
on. This procedure generates an infinite number of fluid equations, which is usually 
truncated at a relatively low order by assuming an appropriate closure equation. 
The closure often takes the form of a thermodynamic equation of state, relating the 
pressure to the density, e.g. the polytropic relation P oc n 1 . The same procedure 
can be applied to the Wigner equation, although some steps in the derivation are 
somewhat subtler than in the classical case. With this technique, quantum fluid 
equations were derived in |29| . Here, we shall present a succinct derivation of 
the same fluid model using a different approach based on the Hartree equations 

[5U1I3T] . 

The starting point is the system of 2N equations Ij4.19p - I|4.20|l . which is com- 
pletely equivalent to Hartree's model (|4 . 1 2J1 - 14.1 3(1 . Let us define the global density 
n(x, t) 

N 

n(x,t) = ^p a n a (4-25) 



and the global average velocity u(x, t) 



N 

u(x,t) = (u a ) = y^P a — u a (4.26) 

a=l 
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By multiplying the continuity equation l|4.19|) by p a and summing over a — 1 . . . N, 
we obtain 

Similarly, for the equation of momentum conservation l|4.2U[l . one obtains 

— +U— = — ^t + — — f p f d xV^ \ 1 dP / 42 g) 

dt dx m dx 2m 2 dx ^— ' \ / mn 9x 

a=l v v 

where the pressure P(x,t) is denned as 

,2 / V m „. \ 2 



P = mn 



= mn{(u 2 a ) - (u a ) 2 ) (4.29) 



So far the derivation is exact, but (|4.28l) still involves a sum over the N states, 
so no simplification was achieved. Our purpose is to obtain a closed system of two 
equations for the global averaged quantities n and u. In order to close the system, 
two approximations are needed: 

1. We postulate a classical equation of state, relating the pressure to the den- 
sity: P = P(n). 

2. We assume that the following substitution is allowed: 

X 'dl.fn 



It can be shown that the second hypothesis is satisfied for length scales larger than 
\p. This will be apparent from the linear theory detailed in Sec. [5j 

With these assumptions, we obtain the following reduced system of fluid equa- 
tions for the global quantities n and u 

ai + -kr = ° (431) 

8 " ^4?^^^)-— (4-32) 



dt dx m dx 2m 2 dx \ y/n J mn dx 

where <p is given by Poisson's equation. We stress that we have transformed a 
system of 2N equations l|4.19l) - (|4.20|) into a system of just two equations. 

An interesting form of the system H4.31|) ~ H4.32|I can be obtained by introducing 
the following 'effective' wave function 

t) = yjn(x, t) exp (iS(x, t)/h), (4.33) 

where S(x,t) is defined according to the relation mu = d x S, and n — \^\ 2 . It is 
easy to show that (|4.31[1 - I|4.32JI is equivalent to the following nonlinear Schrodinger 
equation 

_<9* _ n 2 a 2 * 

d t 2m dx 2 
Wcfr(n) is an effective potential related to the pressure P(n 



^ = -^^-^ + W cS (\m 2 ) *• (4-34) 



i dn' dPM) , . 

W cS (n)=l (4.35) 
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As an example, let us consider a one-dimensional (ID) zero-temperature fermion 
gas, whose equation of state is a polytropic with exponent 7 = 3 



P = -4^ 3 , (4.36) 
3ng 

where vf is the Fermi velocity computed with the equilibrium density no. Using 
the effective potential becomes 



W eS = ^\*\*. (4.37) 



This fluid model is a useful approximation, as it reduces dramatically the com- 
plexity of the Hartree system (2JV equations) or the Wigner equation (phase space 
dynamics). Its validity is limited to systems that are large compared to Xf- Like 
all fluid approximations, it neglects typical kinetic phenomena originating from the 
details of the phase space distribution function. In particular, it cannot reproduce 
Landau's collisionless damping. 

4.4 Vlasov-Poisson. The Wigner and Hartree approaches are both kinetic 
and quantum. The above fluid model was derived by dropping kinetic effects, 
while preserving quantum effects through the Bohm potential. Another possible 
approximation could consist in neglecting quantum effects while keeping kinetic 
ones. The resulting semiclassical limit is, of course, given by the Vlasov equation 

21 + v 2l- -^-^1 - (4 38) 

dt dx m dx dv 

coupled, as usual, to Poisson's equation. 

The Vlasov-Poisson (VP) system has been used to study the dynamics of elec- 
trons in metal clusters and thin metal films [fl I32| . It is appropriate for large ex- 
citation energies, for which the electrons' de Broglie wavelength is relatively small, 
thus reducing the importance of quantum effects in the electron dynamics. How- 
ever, as we have seen in Sec. 12.21 the equilibrium distribution for electrons in metals 
lies deeply in the quantum region (<?q > 1), so that the initial condition must be 
given by a quantum Fermi-Dirac distribution. In this sense, the VP system is 
semiclassical. 

4.5 Initial and boundary conditions. In order to perform numerical sim- 
ulations, boundary and initial conditions must be specified. Physically, the ini- 
tial condition should represent a situation of thermodynamics equilibrium (ground 
state). For electrons in metals, most analytical results are obtained in the case 
of an infinite system (often referred to as the 'bulk' in solid state physics), which 
can be realized in practice by taking periodic boundary conditions with spatial 
period L (this, of course, introduces a lower bound for the wave numbers, namely 
k = 2n/L). For such an infinite system, the ground state can be easily specified. 
For the WP and VP models, any function of the velocity only fo(v) constitutes a 
stationary state. In particular, for fermions we take the Fermi-Dirac distribution 
(see, however, Sec. 15 .21 for a discussion on the appropriate Fermi-Dirac distribution 
for ID problems): 

/ofo) = const, x 1 + e l {e _ ll y (4.39) 

where f3 e = l/fcsT e and e = mv 2 /2 is the single-particle energy (here, v is the 
modulus of the 3D velocity vector). The chemical potential fi(T) is determined 
so that J fodv — uq, where uq is the equilibrium density; [i becomes equal to the 
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Fermi energy Ep in the limit T — > 0. The same ground state can be defined for the 
Hartree model by choosing the initial wave functions in the form 14.22J1 . The Fermi- 
Dirac distribution is then specified by the probabilities p a = [1 + e^ e ( ea_ W] _1 , with 
e a = mul a /2. 

The bulk approximation is not relevant for metal clusters and other metallic 
nanostructures, which are small isolated objects that exist cither in vacuum or em- 
bedded in a background non-metallic matrix. For open quantum systems (Wigner 
or Hartree), the choice of appropriate boundary conditions is a subtle issue, which 
will not be addressed here |33| . For the semiclassical VP system, one can, for 
instance, use open boundaries for the Vlasov equation (zero incoming flux) and 
Dirichlet conditions for the Poisson equation. 

As to the initial condition, the ground state is easily determined for the VP 
system with open boundaries. The main difference from the bulk equilibrium is 
that, for an open system, the electrostatic energy does not vanish and the equilib- 
rium density no(x) is position dependent. As any function of the total energy is a 
stationary solution of the Vlasov equation, we define the initial state as a Fermi- 
Dirac distribution 1)4. 39|) with e(x, v) = mv 2 /2 — e<j>(x), where the potential <j> is 
not yet known. By plugging this Fermi-Dirac distribution into Poisson's equation, 
one obtains a nonlinear equation that can be solved iteratively to obtain <f>, which 
in turns yields the ground state distribution. 

In contrast, stationary solutions of the Wigner equations are not simply given 
by functions of the energy, so that the above procedure cannot be applied. It is 
easier to compute the ground state in terms of the Hartree wave functions ip a , and 
then compute the corresponding Winger function with 1)4. 3[) . Several methods to 
compute the ground state wave functions in the Hartree formalism are available in 
the literature |34L 1351 136j and will not be discussed here. 

5 Linear theory 

In order to compare the various models described in Sec. 01 we shall go into 
some details of the linear theory for a zero-temperature homogeneous equilibrium 
with periodic boundaries, both in ID and in 3D. More detailed calculations on the 
linear theory of quantum plasmas can be found in 2SJ 1371 138| . 

5.1 Zero-temperature ID Fermi-Dirac equilibrium. In one spatial di- 
mension, the Fermi-Dirac distribution at T — is given by fo(v) = no/2vp if 
M < vp and fo(v) — if |v| > vp. The Fermi velocity in ID is 

nhn 

v F = . 5.1) 

2 m 

This distribution is identical to the so-called 'water-bag' distribution, which has 
been extensively used in classical plasma physics 39 . Using the Wigner linear 
dielectric constant (|4.1UI) and developing the results in powers of k and h, one 
obtains the following dispersion relation (details are given in 29 ) 

u 2 = ul + k 2 v 2 F + ^ + M** + {h\ fc 12 ) . (5.2) 
H 4m z 6m z 

With the Vlasov dielectric constant, the dispersion relation for the same equilibrium 

reads as 

u 2 = col + k 2 vp . (5.3) 

Note that the above Vlasov dispersion relation is exact: no terms of order fc 4 or 
higher exist. We also point out that, for this ID equilibrium, the imaginary part 
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of the dielectric constant vanishes identically and therefore there is no Landau 
damping. 

We now want to compare the above result Q5.2JI . obtained with the 'exact' 
Wigner-Poisson model, with the equivalent result obtained with the fluid model 
developed in Sec. 14.31 Let us consider the fluid equations H4.31fH4.32f> m t ne case 
of a zero-temperature ID fermion gas, for which the pressure is given by (|4.36fl . 
Linearizing around the homogeneous equilibrium n — no, u = <f> = 0, one obtains 
the following dispersion relation (no further approximation has been used) 

u*=4 + k\% + ^. (5.4) 

We see that the fluid 15.4fi and the Wigner l|5.2l) dispersion relations coincide up 
to terms of order k 4 . This confirms our conjecture that the fluid model is a good 
approximation of the WP (or Hartree) model in the limit of long wavelengths. 

5.2 Zero-temperature 3D Fermi-Dirac equilibrium. From a physical 
point of view, it is possible to imagine a real physical system that displays ID 
behavior. This can be realized, for instance, in a thin metal film where the dimen- 
sions parallel to the surfaces are much larger than the thickness of the film. The 
electron dynamics can then be described by a ID infinite slab model depending only 
on the co-ordinate x normal to the film surfaces. Even in such a situation, however, 
the ID Fermi-Dirac distribution discussed in Sec. 15. II is not realistic. Indeed, in the 
ground state at zero temperature, the electrons occupy all the available quantum 
states up to the Fermi surface, defined by | v \ < vf- There is no reason why states 
with v y 7^ and v z ^ should not be available, and indeed they are occupied. 
Therefore, the equilibrium distribution is always three-dimensional, even in a ID 
infinite slab geometry: 

/o» = ~4~Y for M ^ "f> and /o» = for \v\ > v F . (5.5) 

3 nV F 

However, it is still possible to keep the ID geometry for the Vlasov and Poisson's 
equations, provided that one uses, as initial condition, the 3D Fermi-Dirac distri- 
bution l|5.5fl projected on the v x axis: fQ B (v x ) = J j fQ D dv y dv z . This yields (we 
now write v for v x ) 

/o 1D («) = ~ (l - J) for M < vf, and f™(v) = for \v\ > v F . (5.6) 

This approach is not as contrived as it might appear at first sight. Indeed, 
linear wave propagation in a collisionless plasma is intrinsically a ID phenomenon, 
which involves plane waves traveling in a well-defined direction. In computing the 
dispersion relation from the 3D equivalent of l|4.10fi or l|4.11fi . we would first inte- 
grate over the two velocity dimensions normal to the direction of wave propagation 
(which can, arbitrarily and without loss of generality, be chosen to be the v x direc- 
tion). Therefore, we would be left with a reduced distribution such as (|5.6fl that 
intervenes in a ID dispersion relation such as l|4.10f) or (|4.11fi . This line of reason- 
ing is no more valid when nonlinear effects become important, as these may trigger 
truly 3D phenomena. Collisions also constitute an intrinsically 3D effect. 

We now insert 1)5. 6fl into (|4.1Ufi or (|4.11l) and compute the dispersion relations 
for the WP and VP systems, developed in powers of k. For the equilibrium l|5.6fl . 
the dielectric constant does display an imaginary part, and therefore there is, in 
principle, some collisionless damping associated with this equilibrium; however, we 
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shall neglect it for the moment and concentrate on the real part of e(k, ui). We also 
assume the following ordering 

hk LU 

— < v F < -, (5.7) 
to k 

which is valid for long wavelengths. Note that the second inequality in 15.7(1 means 
that the phase velocity of the wave must be greater than the Fermi velocity. With 
this assumption, the dispersion relation up to fourth order in k reads as 

J=ul + \ k 2 v% + ^ + . ... (5.8) 

Let us now derive the dispersion relation for the quantum fluid model 14.31JI - 
(|4.32|) . by assuming an equation of state of the form: 

where Pq and no are the equilibrium pressure and density, respectively. We obtain: 

^=a£+ 7 k a «S + ^, (5.10) 

where Vq = Po/(mno). We note that the ID fluid dispersion relation 15. 4|) is 
recovered correctly when 7 = 3 and Pq = tiq mv F /3 as can be deduced from (|4.36H . 
Now, in 3D, the pressure of a quantum electron gas at thermal equilibrium and 
zero temperature can be written as |13| 

Po = ln E F (5.11) 
5 

(where Ep is computed at equilibrium), so that Vq = vp/5 and l|5.10[l becomes 

- 2 =-p + J^o + SJ • (5-12) 

One may think that the correct exponent to use in the equation of state (|5.9() is 
7 = (D + 2)/D, yielding 5/3 in three dimensions. However, with this choice, the 
fluid dispersion relation (|5.12|) would differ from the WP result 15. 8|) . The correct 
result is obtained by taking 7 = 3, just like in the ID case. Why is it so? The reason, 
as was already mentioned earlier, is that linear wave propagation is essentially a ID 
phenomenon, because it involves propagation along a single direction, without any 
energy exchanges in the other two directions. The details of the linear dynamics are 
essentially determined by the equation of state, which must therefore feature the 
ID exponent. In contrast, the equilibrium is truly 3D (because we have projected 
the 3D Fermi-Dirac distribution over the x direction): therefore, the equilibrium 
pressure must indeed be given by its 3D expression (|5.11|) . 

5.3 Finite temperature solutions: Landau damping. So far, we have 
completely neglected the fact, discovered by Landau in 1940 40 ; , that electrostatic 
waves can be exponentially damped even in the absence of any collisions. The rig- 
orous theory of Landau damping can be found in most plasma physics textbooks. 
Here, we shall only remind the reader that the damping originates from the singu- 
larity appearing in the dispersion relations l|4. 9fl and (|4.11|) at the point v = ui/k in 
velocity space. This corresponds to particles whose velocity is equal to the phase 
velocity of the wave tu/k (resonant particles). Landau showed that the correct way 
to perform the integral in the dispersion relation is not simply to take the prin- 
cipal value (which only yields the real part of the frequency) , but to integrate in 
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Figure 3 Fermi-Dirac velocity distribution projected on a single velocity di- 
rection, with T/T F = 0.05. 



the complex v plane, following a contour that leaves the singularity always on the 
same side. With this prescription, the dielectric function is found to possess an 
imaginary part, which in turn gives rise to a damping rate for the wave. This 
argument, originally developed for the VP system, still holds for the quantum WP 
case, although of course the numerical value of the damping rate will depend on h. 

At zero temperature, no particles exist with velocity v > vf- Therefore, waves 
with phase velocity larger than vf are not damped at all. For these waves, we have 
k < uj/vp', as the real part of the frequency is approximately equal to the plasma 
frequency, this means that waves for which k\p < 1 are not damped. These are 
waves with a wavelength smaller than = vf l^pe-, which is of the order of the 
Angstrom for metallic electrons (see table 1). At finite temperature, the projected 
Fermi-Dirac distribution i|5.6[l (see Fig. |3J extends to all velocities (although it 
decays quickly), so that some amount of damping exists for all wave numbers. 

The linear damping rate can be computed from the dispersion relation, see for 
instance |41| , and we shall not deal with this issue any further. It is more interesting 
to look for some qualitative guess about the nonlinear phase that follows the initial 
Landau damping |42| . Classically, Landau damping lasts up to times of the order of 
the so-called bounce time m,, after which the evolution is inherently nonlinear. The 
bounce time is related to the amplitude of the initial perturbation of the equilibrium 
distribution 

f(x,v,t = 0+) = f (v) (1 + a cos kx), (5.13) 

where a = n/ny is the normalized density perturbation and k is the wave number 
of the perturbation. The bounce time can then be written as: ujpTt, — or 1 / 2 . 
When the perturbation amplitude is not too small, one generally observes that 
Landau damping stops after a time of the order r\,. This happens because resonant 
particles (whose velocity is close to the phase velocity of the wave) get trapped 
inside the travelling wave, thus creating self-sustaining vortices in the phase space. 
The presence of such vortices maintains the electric field to a finite (albeit generally 
small) level. 
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Figure 4 Classical phase space portrait of the electron distribution function 
in the region around the phase velocity of the wave. Position is normalized to 
\p and velocity to vp. The simulation parameters are: T = O.OlTp, a = 0.1, 
and k\p = 1. 



We have performed numerical simulations of the VP system using a Vlasov 
Eulerian code |43l I44| . The initial equilibrium is given by the projected Fermi- 
Dirac distribution at finite temperature 



ft ^ 3n ° T i 
feO(X,V x ) = - — —In 

4 vp lp 



1 + exp 



(5.14) 



which generalizes the zero-temperature result (I5.6|) . We took an equilibrium tem- 
perature T = 0.01 Tp, whereas the perturbation l|5.13[) is characterized by an 
amplitude a = 0.1 and a wave number k\p = 1. The phase-space portraits of the 
electron distribution (Fig. 0} show, as expected, the formation of a vortex traveling 
with a velocity close to co/k. Further, it can be easily proven that the extension of 
the vortex in velocity space Sv is related to the perturbation's amplitude and wave 
number in the following way 



5t, = ^ a Va. (5.15) 
k 

We now turn to the fully quantum case, described by the Wigner-Poisson sys- 
tem with the same initial condition (|5.13l) . The simulations have been performed 
with the code described in Ref. |41| . As the wavelength of the perturbation is 
2-7r/fc, the classical phase-space vortex defines a phase space area of order Sv/k. If 
this area is smaller than H/m, then quantum mechanics forbids the creation and 
persistence of such a structure 1411 . Using the relation l|5.15|l . we find therefore 
that the phase-space vortex should be suppressed by quantum effects when 

— > ^ a 1/2 . (5.16) 

m k 

Using normalized units, the above relation becomes 

H k 2 X 2 F > v^, (5.17) 

where H = huj p /Ep is a measure of the magnitude of quantum effects. As quantum 
effects prevent particles from being trapped inside the wave, we expect Landau 
damping to continue even for times larger than the bounce time. 
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Figure 5 Quantum phase space portrait of the electron Wigner function in 
the region around the phase velocity of the wave. Position is normalized to 
\p and velocity to vp. Same parameters as in Fig. [3] with in addition H = 1. 
Left frame: positive part of f(x,v); right frame: negative part of /. 




Physically, the above effect is related to quantum tunnelling: particles trapped 
in the wave have a certain probability to be de-trapped, even if their energy is less 
than that of the potential well of the wave. Yet another way to picture this effect is 
to consider that, if the potential well is too shallow (i.e. for small a), no quantum 
bound states can exist inside it. 

We have tested the above order-of-magnitude arguments by running the same 
simulation as that of Fig. 0] but using the Wigner, instead of Vlasov, equation. 
We take H = 1, so that the inequality (|5.17(l is satisfied. The phase space portraits 
(Fig. [SJ indeed confirm that no vortex structures appear. Note also the appearance 
of large areas of phase space where the Wigner distribution function is negative. 
Comparing the classical and quantum evolution of the potential energy (Fig. EJ, 
we observe that, for long times, the electric field is significantly smaller in the 
quantum evolution. These results suggest that semiclassical models of the dynamics 
of metallic electrons may underestimate the importance of Landau damping. 
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6 Conclusions and future developments 

In this paper, we have reviewed a number of kinetic and fluid models that are 
appropriate to study the behavior of collisionless plasmas when quantum effects 
are not negligible. The main application of these models concerns the dynamics of 
electrons in ordinary metals as well as metallic nanostructures (clusters, nanopar- 
ticles, thin films). In particular, it is now possible to study experimentally the 
electron dynamics on ultrafast (femtosecond) time scales that correspond to the 
typical collective plasma effects in the electron gas. The properties of a quantum 
electron plasma (neutralized by the background ions) can thus be measured with 
good accuracy and compared to the theoretical predictions. 

The main drawback of the models described in this paper is that they all ne- 
glect electron-electron collisions. Strictly speaking, the collisionless approximation 
should be valid only when gQ -C 1, which is not true for electrons in metals (see 
Table 1). As we discussed in Sec. Pauli blocking does reduce the effect of col- 
lisions for distributions near equilibrium, but many interesting phenomena involve 
nonequilibrium electrons, so that the validity of the collisionless approximation is 
not completely clear. It is possible, in principle, to include collisional effects in 
semiclassical models such as the Vlasov-Poisson system, by simply adding a col- 
lision operator on the right-hand side of the Vlasov equation p]. The collision 
operator relevant for a degenerate fermion gas is known as the Uhling-Uhlenbeck 
collision term and is basically a Boltzmann collision operator that takes into ac- 
count the exclusion principle. Even in the semiclassical case, however, the validity 
of the Uhling-Uhlenbeck approach may be questioned for strongly coupled plasmas 
(gQ > 1), for which it is conceptually difficult to separate the mean-field from the 
collisional effects. 

It is much harder to include a simple collisional term in truly quantum models, 
such as the Wigner or Hartree equations [45] . There is a vast literature on dissi- 
pative quantum mechanics, but this is concerned with the interaction of a single 
quantum particle with an external environment |46l 1471 148) . In addition, the cou- 
pling to the environment is generally assumed to be weak, which is not the case for 
metallic electrons. Recent approaches to the dynamics of strongly coupled quantum 
plasmas range from quantum Monte Carlo methods (for the equilibrium) to quan- 
tum molecular dynamics simulations (for the nonequilibrium dynamics) |15L 149] . 

We shall mention another two possible extensions of the models presented in 
this paper, but these are more technical points, rather than conceptual issues such as 
the above-mentioned inclusion of electron-electron collisions in the strongly coupled 
regime. 

The first issue concerns the coupling to the ion lattice. So far, we have as- 
sumed that the ions form a motionless positively-charged background described by 
the equilibrium density no(x), which may be position-dependent for open systems 
such as clusters. This is appropriate for times shorter than the typical ion-electron 
collision time (see Table 1), but for longer times the ion dynamics must be taken into 
account. The ion dynamics, just like the electrons', may be split into a mean-field 
part and a collisional part. If we want to include only the mean-field component 
(and neglect ion-electron collisions), then all we need to do is add a Vlasov equation 
for the ionic species, supplemented by a Maxwell-Boltzmann initial condition (as 
the ions are always classical) - see, for instance, |32j . This is conceptually simple, 
though it may require rather long simulation times, because the typical ionic time 
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scale (proportional to u> P i) is much longer than the electrons'. As a first approxi- 
mation, electron-ion collisions may be modelled by a simple relaxation term of the 
Bathnagar-Gross-Krook type 50 , which is the kinetic analog of the Drude relax- 
ation model of solid state physics |14| . For a more accurate approach that includes 
the full ion dynamics (mean-field and collisional), it will be necessary to treat the 
ions with molecular dynamics simulations. 

Finally, one should consider the effect of magnetic fields (both external and 
self-consistent) on the plasma dynamics. Magnetic fields should not alter the main 
conclusions drawn in the present work and are easily included in all the equations 
presented here. For very strong fields, fusion and space plasma physicists have 
developed a battery of approximations (guiding-center, gyrokinetic, . . . ) that allow 
one to reduce the complexity of the relevant models. It is a challenging task to 
transpose these methods to the physics of quantum plasmas |51L 152) . Magnetic 
fields should also trigger spin effects, which are uniquely quantum-mechanical. The 
interaction of spin and Coulomb effects, still a largely unexplored field, is thus likely 
to become an active area of future research. 
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